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Abstract 

We present a manifestly covariant quantization procedure based on the de Donder-Weyl Hamil- 
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CN \ tonian formulation of classical field theory. This procedure agrees with conventional canonical 

el ■ 

quantization only if the parameter space is d = 1 dimensional time. In d > 1 quantization requires 
. a fundamental length scale, and any bosonic field generates a spinorial wave function, leading to 

{Sj ' the purely quantum-theoretical emergence of spinors as a byproduct. We provide a probabilistic 

>' 

Q^ ' interpretation of the wave functions for the fields, and apply the formalism to a number of simple 

on : 

examples. These show that covariant canonical quantization produces both the Klein-Gordon and 

0^ ; 

the Dirac equation, while also predicting the existence of discrete towers of identically charged 
fermions with different masses. Covariant canonical quantization can thus be understood as a 

^ : 

• 'first' or pre-quantization within the framework of conventional QFT. 
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I. INTRODUCTION 



The apparent incompatibility between general relativity and quantum mechanics has long 
been a topic of concern and interest in the theoretical physics community. Diffeomorphism 
invariance has to be satisfied on the side of a general relativistic theory, in particular denying 
any fundamental distinction between the notions of space and time; but it is less clear how to 
achieve this requirement in, or properly translate it to, a quantum theory. This particularly 
applies to the canonical formulation of quantum mechanics and quantum field theory based 
on a Hamiltonian treatment. A neat way around this problem may be seen in path integral 
quantization which explains why the predictions of the quantized theory still possess the 
relativistic symmetries of the classical theory; but from the Hamiltonian point of view with 
its explicit space-time split this is not a special merit of the quantization procedure. This 
motivates the question whether there is a covariant extension of Hamiltonian methods which 
also allows for a manifestly covariant quantization procedure. 

On the level of classical field theory there is indeed a Hamiltonian formulation that does 
not rely on singling out a time coordinate, but treats all spacetime coordinates equally 
throughout. This theory was presented already in the nineteen-thirties by de Donder l| 



and Weyl |2j. Full covariance is maintained through the use of multi-momenta, where one 
momentum is associated to each partial derivative of the fields. While providing a fully 
covariant equivalent to the standard Hamiltonian formulation of field theory (in the sense 
of providing the same solutions), the de Donder- Weyl formulation of classical dynamics 
has not received too much attention. Only recently have there been several attempts to 
quantize field theories on its basis. An early attempt by Good has been shown to 

disagree with ordinary quantum mechanics and to give incorrect predictions for the hydrogen 
spectrum p. Subsequently, a quantum equation based on de Donder- Weyl theory has been 



conjectured by Kanatchikov 6 1 and Navarro [7J. There have also been attem pts to obtain 
a path integral formulation [g, aQ] and a version of Bohmian mechanics jlll I12I based on 
de Donder- Weyl dynamics. Other recent applications of the de Donder-Weyl formulation 
of field theory include a derivation of the Ashtekar-Wheeler-DeWitt equation of canonical 
quantum gravity 

In this paper, after a brief review of some of the elements of the classical de Donder- 
Weyl theory in the following section we will formulate a covariant Poisson bracket. In 
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section HTT1 we will then proceed to apply the Dirac quantization postulate to the latter. If 
supplemented with a second, geometrically motivated, quantization postulate, this leads to 
the same quantum evolution equation that had previously been conjectured on the basis 
of analogies [3, Q|. Our approach for the first time presents a derivation of this equation, 
which unifies both the Schrodinger and the Dirac equation, from first principles. We go on 
to develop the quantum theory in the covariant Schrodinger picture; in particular, we will 
discuss the representation of operators, the consequences of an indefinite scalar product on 
the Hilbert space which immediately follows from the requirements of covariance, and the 
probability interpretation of the wave functions. We apply the theory to a number of basic 
problems in section \W\ with sometimes surprising results: 

Among them are a new derivation of the Klein-Gordon equation that makes no use of the 
relativistic energy-momentum relation, the emergence of spinors from the quantization of 
scalar theories, and in particular the emergence of the Dirac equation from the quantization 
of any scalar field action. This means that the quantization procedure here presented does 
not replace quantum field theory; instead, it is found to provide a supplementary "first 
quantization". A result of potential phenomenological interest is the prediction of towers 
of identically charged fermions that differ only by their masses, providing a qualitative 
explanation for the generations in the Standard Model. 



II. COVARIANT HAMILTON! ANS IN CLASSICAL FIELD THEORY 

This section reviews the covariant Hamiltonian treatment of classical field theories, dis- 
cussed first by de Donder and Weyl , which is based on the introduction of multi- 
momenta associated to the partial derivatives of the fields. We then define a new covariant 
Poisson bracket to rewrite the general phase space evolution equations in an equivalent form 
suitable for quantization. 



A. Field theory in the multi-symplectic formalism 

Consider a geometrically well-defined field theory which is diffeomorphism invariant on 
a (i-dimensional Lorentzian background manifold E. We will call this background manifold 
the parameter space of the theory, coordinatized by parameters {er a } with corresponding 
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partial derivatives d a = d/da a . Classical fields q % are functions on this manifold, i.e., 

q* : E -> R. (1) 

In some theories it is convenient to consider a set of n fields {g 4 } as coordinates of a second, n- 
dimensional, target space manifold M; in this case, requiring the theory to be geometrically 
well-defined means it should obey the further diffeomorphism invariance on M. The notion 
of a target space manifold is, however, secondary. We define the theory on E by its action, 
which is obtained from the integration over E of a scalar Lagrangian as 

S = J^d d o^L(q\V a qi). (2) 

Note that the standard quadratic kinetic term in the Lagrangian depends on V a q l - Forming 
a scalar from these (covariant) derivatives necessitates the existence of a non-degenerate, 
and hence invertible, metric g on E, the signature of which we take to be (— , +,..., +). The 
determinant of this metric appears in the integration measure. An explicit dependence of L 
on the coordinates of E is excluded by the requirement of diffeomorphism invariance. 

The equations of motion of the theory (J2J) are the Euler-Lagrange equations derived by 
variation of the action with respect to the fields, 

" 81 £=0, (3) 



a dV a q i dq i 

where the covariant derivative involves the unique torsion free and metric compatible Levi- 
Civita connection of g. The necessary boundary condition requires a vanishing integral 

L dS W^°- (4) 

The above equations of motion are partial differential equations of second order, for first 
order Lagrangians. To reduce the order, the standard Hamiltonian treatment introduces 
canonical momenta pi = dL/ddoq 1 . Clearly, these momenta are non-covariant quantities as 
their definition explicitly depends on the choice of time and hence on the choice of coordinate 
system on E. It follows that the usual Hamiltonian function, depending on the non-covariant 
canonical momenta, cannot be a scalar. 

To remedy this apparent difficulty, we introduce the manifestly covariant multi-momenta 
associated to each partial derivative of the fields, 



which transform as the components of a vector in the parameter space tangent bundle TE 
(and as those of a differential form in T*M, if the fields form coordinates of a target space 
manifold). We assume Lagrangians such that the multi-momenta as functions of the fields 
and their partial derivatives may be solved for these derivatives to yield V ' a <t (q^pl)- I n 
terms of the new covariant momenta, we may then also define the covariant Hamiltonian 

H = P y a q l ~ L (6) 

which is a function of the new independent variables q % and p", and transforms as a diffeo- 
morphism scalar on the parameter space E. 

The Euler-Lagrange equations imply the covariant Hamiltonian equations 

dH 

^ = -VaP? , (7a) 
f) M 

Conversely, given a covariant Hamiltonian H (q\p^), we may define a Lagrangian 
L(g*,V a g : ') via ©. Then the covariant Hamiltonian equations imply the Euler-Lagrange 
equations. Diffeomorphism invariance again implies that the Hamiltonian cannot depend 
explicitly on the coordinates of S. Below we will see that the covariant Hamiltonian for- 
malism nicely reduces to conventional Hamiltonian mechanics if the parameter space E is 
one-dimensional . 



B. The classical Dirac field as an example 

As an example for the powerful finite-dimensional phase space formalism of de Donder 
and Weyl, we take a brief look at the massive Dirac field. The Lagrangian in its symmetrical 
form is given by 

L = l*Pl a V a *P - ~ V^tV - Mipip , (8) 

where we have introduced the Dirac matrices 7° of the (curved) background, on which we will 
comment in more detail below. We treat i/j and ip as independent, such that the conjugate 
covariant momenta follow by definition as 

K = \^ a , vr| = -i 7 >. (9) 
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These relations are, in fact, primary constraints, relating the spinors and their conjugate 
momenta. Although these momenta are not invertible to obtain V a ip(ip,tp, 7T^,7Tt), and 
similarly Va0, we can define the covariant Hamiltonian as 

H = + (vr; - fy a /2) A a + A a (tt| + 7^/2) , (10) 

where the constraints have been added with the help of spinorial Lagrange multipliers A a 
and A a . The Dirac equation and its conjugate now follow immediately from the covariant 
Hamiltonian equations (JJJ) above, utilizing the constraint equations. 



C. Definition of a covariant Poisson bracket 



With the aim of facilitating an easier transition to a quantum theory, we consider Poisson 
brackets in the new formalism. A covariant extension of the standard Poisson bracket is given 
by the definition 

If 1 - ^L^l_ _ d£dg_ 

V>9la Q q iQ p a Q p a Q q i \ > 

for any two phase space functions / and g depending on the Hamiltonian variables q l and p". 
This bracket carries a further index, thus mapping two functions of the canonical variables 
to a differential form in T*E. In general, it changes the number of indices and with it the 
tensor structure defined by its arguments. This obstructs the usefulness of this bracket 
definition, as valuable properties of the Poisson bracket are lost. This applies in particular 
to the important Jacobi identity, which provides the algebra of phase space functions with 
the structure of a Lie algebra. Here the Jacobi identity is valid only for equal subscripts, 
i.e., for expressions of the form {{/, g} a , h} a , but these are not allowed as tensors on E. 
Hence we are led to amending the bracket definition, and consider brackets of the form 

{f,g} = {f,g} a t a (12) 

where we introduce an arbitrary vector field in TE with components t a , the origin of which 
we will discuss in the following section on quantization. The only classical requirement that 
we will make on this field concerns its normalization N(d) = g a bt a t b which may depend on 
the dimension of the parameter spacetime E. It should be such that N(l) = — 1 so that the 
usual Poisson bracket may emerge when d = 1 with g aa = — 1. It turns out that the bracket 
so defined satisfies the formal algebraic properties of the Poisson bracket, which we state for 



phase space functions / and g and real numbers c[ 



within the de Donder-Weyl formalism, compare 
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or other definitions of Poisson brackets 



16 



m 



The Poisson bracket is antisymmetric and annihilates constants, 

{f,g} = -{g,f}, (i3a) 

{/,c} = 0; (13b) 
it is IR-linear in / and g (where linearity in the second argument follows by antisymmetry), 

{h + h,g} = {h,g} + {f 2 ,g}, (14a) 

{cf,g} = c{f,g}; (14b) 
the Poisson bracket further satisfies a product rule and, importantly, the Jacobi identity: 

{fif2,g} = {fi,g}f 2 + A {/*</}, (15a) 
{{/, g},h} + {{g, h},f} + {{h, f},g} = Q. (15b) 

The introduction of the vector field t a provides another advantage, again with a view 
towards quantization: it allows us to achieve a one to one correspondence between the 
fields q % and the contracted multi-momenta Pi = —t a p^. We find the covariant Poisson 
brackets 

{q\q 3 } = 0, (16a) 
{Pi,Pj} = 0, (16b) 
{q\p 3 } = -N(d)S]. (16c) 

Brackets including the covariant Hamiltonian generate the following expressions, similar 
to those appearing in the covariant Hamiltonian equations of motion ((7j); the fact that there 
is no precise agreement is due to the appearance of the vector field t a : 

{q\H} = t a V a q l , (17a) 
{ Pl ,H} = -N{d)V a pt. (17b) 

These brackets are useful in evaluating the evolution equation for phase space functions with 
respect to the parameters given by the coordinates of E. We calculate 



which may be rewritten as 



ev a f-{f,H}-t a v°j 



df_ 

dp b i 



- h (t a V aP \-{plH}). 



(19) 



This is the form of the general evolution equation which we will use as an important ingre- 
dient of the quantization procedure. Note that the parameter space derivatives of the phase 
space function / are evaluated along the integral curves of the vector field t a . The derivative 
operator V° acts only on the cr a -dependence of / not coming in through the coordinates 
and momenta. A closer inspection of the equation also reveals that it is trivially satisfied 
for any phase space function linear in the momenta, e.g., for pi or t a pf. This means we have 
to supplement it with equation (|17b[) . 

D. Hamiltonian mechanics on one-dimensional £ 

The results and constructions above are in complete analogy to the standard Hamiltonian 
treatment of classical mechanics which is, however, restricted to a one-dimensional parameter 
space E, with time coordinate a, if diffeomorphism invariance is required. 

The Hamiltonian formalism of de Donder and Weyl reduces to the standard one for d = 1. 
To see this more explicitly, note that the TE index a of the multi-momentum pf can merely 
take a single value in this case corresponding to the single coordinate o on E, which may be 
suppressed. The manifold, its tangent and cotangent spaces are all locally isomorphic to the 
real numbers. The normalization requirement for the single-component vector field enforces 
t u = 1 because of our signature convention g aa = — 1. Thus we obtain agreement between 
our covariant Poisson bracket and the standard one. Equations ()16|) reduce to the canonical 
Poisson brackets, and equations (fTTjl become equivalent to the Hamiltonian equations. The 
right hand side of the phase space evolution equation (fiT)|) cancels; what remains is the 
well-known time evolution formula 



d_ 

da 



f-{f,H}-— f = 0. 



(20) 



III. COVARIANT CANONICAL QUANTIZATION 



Following Dirac, the quantization of a system in classical mechanics takes as its starting 
point the Hamiltonian formulation. The canonical variables are promoted to operators 
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acting on a Hilbert space, and the Poisson brackets to commutators. With our new covariant 
Hamiltonian formalism, we will now mimic these steps. 

A. Quantization postulates 

As the bracket {/,<?} that we have defined above has the same algebraic properties as 
the Poisson bracket, we promote it to a commutator of operators in exactly the same way, 

{f,g} ^ -U^ilg}- (21) 

According to a famous argument of Dirac jlfij . this is in fact the only consistent quantization 
postulate, if the quantum bracket is required to preserve its classical algebraic properties. 
Phase space functions / and g have been replaced by operators on some Hilbert space, 
denoted by a hat. The imaginary unit is required to imply that — i [•, •] is self-adjoint for self- 
adjoint entries (with respect to the Hilbert space inner product which we will define below). 
In our units, where c = 1 and h = 1, we have to introduce another independent length scale 
/, which we might choose to be the Planck length lp, to compensate the dimensions of the 
derivatives with respect to the canonical variables that appear in our Poisson bracket. This 
result does not depend on the dimension of the fields q 1 ; it merely assumes the dimension 
of the Lagrangian L is (length) ~ d . Note that the necessity of a fundamental length scale for 
quantization appears only on manifolds £ of dimension d > 1. 

Now we have to think about the vector field t a in our bracket definition. The obvious 
choice seems to be a classical timelike vector field on E, with normalization N(d) = — 1 
for any dimension. However, this would have several undesirable consequences; firstly, no 
such vector field was included in the classical theory in the original formulation (J2J). Thus 
additional input would be necessary for the quantum theory. Such input would not be 
universal in the sense that the chosen timelike field could differ for different quantum systems 
under consideration. Secondly, this would amount to introducing a space-time split of £ into 
the product of a family of timelike curves and their corresponding normal surfaces, thereby 
introducing all of the problems associated with the canonical procedure. Thus the aim of 
an intrinsically higher-dimensional quantization procedure on £ would be lost. But what 
options are left now of choosing a vector field which is implicitly given on any E? 

To answer this question, we have to resort to some more geometry. On every curved 
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background manifold X (admitting a spin structure) exists an algebra of Dirac matrices 7 a 
with the property that 

lalb + 767a = 2 #a6 • (22) 

As usual, these Dirac matrices are related to those of the local Lorentzian tangent spaces 
by the vielbeins as 7 a = e^T^. The vielbeins form the metric as g a b = e^e^T]^. The 
normalization of the Dirac matrices gives 7 a 7 a = d. This leads us to propose the following 
quantization postulate for the vector field t a , 

t a ^ -m(rf) 7 a , (23) 

which implies N(d) = —dn(d) 2 . To satisfy the requirements N(l) = —1 and t u = 1 for the 
one-dimensional case, we note that the only Dirac matrix in d — 1 is P 7 = i, so that the 
normalization function n(d) must be chosen such that n(l) = 1. Otherwise n(d) is quite 
arbitrary and must be fixed by application of the theory, which remark also applies to the 
fundamental length scale I. 

It is worth noting that the appearance of the Dirac matrices in this context has a historical 



parallel in Dirac's original derivation of the Dirac equation 
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, where a universal object 



with a covariant vector index was also required to fulfill the demands of covariance. 



B. Quantum evolution — Dirac is Schrodinger 

We will now motivate a quantum evolution equation based on our classical covariant 
Hamiltonian picture, which turns out to unify the Dirac and the Schrodinger equation. 

The multi-symplectic phase space is spanned by the canonical variables q l and pf. Any 
'proper' phase space function / depends on the coordinates of the parameter space S only 
through these variables. For the operators / associated to such phase space functions we 
now analyze the requirement that the classical evolution equation (|T9"j) holds in its quantum 
version as 

l l d - 1 [f,H] -in(d) (y /) =S\J^_ (il^iplH] -in(d) {ffflfj . (24) 

The operator S denotes a symmetrization of the canonical variables, which are now opera- 
tors, and we use Feynman's shorthand notation y = 7 a V a . Noting that (y/) = [y , /] in 
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the action on states, the above equation is can be rewritten in the form 

[f,H + n(d)r d+1 f] =S (^-[piH + n {d)l- d+1 y\^ . (25) 

This equation holds for all operators if, and only if, H+n(d)l~ d+l y is a constant independent 
of the canonical variables. For d = 1, this reduces to H + id a . To obtain the same result as in 
conventional quantum mechanics in this limit, we have to set H +n{d)l^ d+1 ^7 = 0. However, 
note that this type of derivation in quantum mechanics does not produce the Schrodinger 
equation, which is H — id a = 0, acting on Schrodinger picture states. This is because the 
quantization of the classical evolution equation leads to an operator equation valid in the 
Heisenberg picture. The unitary change of pictures is responsible for the change of signs. 

We assume, as will be justified in the following section UlI C\ as it requires some develop- 
ment of the theory, that the same change of signs occurs here. Thus we finally arrive at the 
quantum evolution equation, as an equation acting on ^-dependent states: 

(H - n(d)r d+1 f ) |^((T a )) = . (26) 

Note that the quantum evolution effectively does not require the terms appearing in the 
symmetrization operator, which thus need not be specified. This is an advantage because 
it is not consistently possible to do so even in conventional quantum mechanics: no map of 
phase space functions into an operator algebra exists, compatible with the Poisson bracket. 
An attempt to rectify this situation is made by deformation qua ntization, employing as 
operators formal power series in the Planck quantum h, see 21, 22] for reviews. 

The above equation gives the operator H the dimension of mass times I +1 . Supposing, 
in an expansion in terms of the canonical variables, that there is a constant term in H, we 
find that both the Schrodinger and the Dirac equation follow from the same quantization 
procedure. The Schrodinger equation is relevant for a quantization of fields x l : £ = K. — > M n 
and the Dirac equation corresponds to quantized fields on a parameter spacetime £ = M\^. 
This will be further illustrated below. Here we only note that, in d > 1, the wave functions 
will automatically become spinors. 

The quantum evolution equation in the form ()26|) has been conjectured before by 
Kanatchikov and Navarro [2^ on the basis of analogies between the Dirac equation 



and conventional quantum mechanics. Here we have presented for the first time a physi- 
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cally reasonable derivation of equation (|2l)j) from first principles, based solely on the two 
quantization postulates (|2Tj) and 



C. The local evolution operator and the Heisenberg picture 

We shall now introduce the covariant Heisenberg picture, and in this context we will 
justify the remaining assumptions going into the derivation of the quantum evolution equa- 
tion ()26p. To clarify the calculations we will use subscripts g and # for Schrodinger and 
Heisenberg picture quantities, respectively. 

The first notion we need is that of an evolution operator U(a,ao), with the help of 
which a E-dependent Schrodinger picture state may be written in terms of a E-independent 
Heisenberg picture state: 

\^{a)) s = U{a,a Q )\^{a Q )) H . (27) 

The consistency of expectation values requires U to be unitary in the sense U#U = 1 (cf. 
section IIII El below) . 

On curved parameter spaces E, such an evolution operator can only be defined locally, 
i.e., in a sufficiently small neighborhood of a point p e E. If q is another point in this 
neighborhood, then there is a unique geodesic joining p and q. The corresponding evolu- 
tion operator U(q,p) can be written as U(a,ao) i n a geodesic normal coordinate system. 
Substituting the above state expansion into the quantum evolution equation we find 

H s U(a, a Q ) - n{d)r d+1 y U(a, a ) = . (28) 

Solving this equation to first order in an infinitesimally small displacement 5a then gives 

U(a + 5a, a ) = l + ^—5a a la H s . (29) 

The adjoint £/ # , again to first order, follows from (j a Hs)^ = —Hsja- Hence U is unitary, 
if, and only, if 

[7 a ,^s] = 0, (30) 

which we must require. (The examples below show that this relation usually holds.) The 
evolution operator for finite coordinate differences within the local neighborhood follows 
from the limiting procedure U(a, a ) = lim (U(a, a — 5a) . . . U(a + 5a, a )) for 5a — > 0; in 
consequence, it is unitary as well. 
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We are now in the position to calculate the equation of motion for Heisenberg picture 
operators /# = U#fsU- Using the same steps as in conventional quantum mechanics yields 



^/ f H -n{d)- 1 l d - l [f Hl H H ] = [V a f> # ,7 a ]t>/H. (31) 

The right hand side of this equation does not contribute. This is easily seen by noting 
that [V a £7 # ,7 a ] = [V tt f7, 7 a ] # , and using equation to determine V a U ~ l a Hs- Thus 
one finds [V a U, 7 a ] ~ [la,Hs], the vanishing of which was required by the existence of the 
Heisenberg-Picture. The Heisenberg equation of motion 

n(d)r d+1 ff H -[f H ,H H ]=0 (32) 

follows. Using the same method as in the derivation of the quantum evolution equation in 
section UlI Bl one thus finds Hh + n(d)l~ d+1 y = in the Heisenberg picture, fully justifying 
the change of sign we made in obtaining equation (|2fij) . 

The Heisenberg equation of motion (|3*2*|) will play an important role in showing that our 
theory indeed has the correct classical limit in section UlI Fl below. 



D. Canonical operators in the Schrodinger picture 

In the covariant Schrodinger picture, all operators, including the covariant Hamiltonian, 
act on states which are elements of some Hilbert space TC and depend on the coordinates of 
S, i.e. \ip(a a )). We wish to find a realization of these operators acting on wave functions in 
an explicit Schrodinger representation; for this purpose we have to introduce a basis of Tl, 
which is conveniently given by states 

\q a ) = \q) <8> e a , (33) 

where \q) are eigenstates of the field operators q\ i.e. q l \q) = q l \q), and e Q are the canonical 
basis vectors of the representation space of the Dirac algebra. The dual states are given by 

( a q\ = (q\®ui a ) (34) 

where {uj 01 } is the dual basis of {e a }, such that the normalization condition becomes 

( a q | qp) = 6%5(q - q) . (35) 
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In the basis any state of the Hilbert space can be expanded as 

dq\q a )r(q)- (36) 



The components ip a {q) = ( a q \ ip) with respect to this basis give the spinorial Schrodinger 
picture wave function. In this notation we have suppressed the ^-dependence; more precisely, 
one should write ip a {°' a ] <?)• The identity operator on 7i has a partition of the form 



1 = J dq\q a )( a q\ . (37) 
The canonical operators should satisfy the commutation relations 

W,q j } = 0, (38a) 
\pi,Pj] = 0, (38b) 
= idn{d)H- d+l 5), (38c) 

which follow from an application of the two quantization postulates to the classical Pois- 
son bracket equations ([16)1 . It would seem to be convenient at this stage to remove the 
dimension dependence of the canonical commutation relations by setting n(d) = 1/y/d, but 
we emphasize again that n(d) should be fixed by application. To study the action of the 
canonical operators on wave functions, we need the following matrix elements: 



d_ 

( a q I l a I Op) = (lT^(q - q) , (39b) 



*q\Pi\Qf>) = idn(d) 2 l- d+l 5 a p -^-5(q-q), (39a,) 



where the first identity follows from an expansion of ( a q \ [q\pi] \ qp), using ()38cj) . Now we 
act with our operators on arbitrary states, which yields 



= J dq\q a ) (qr(q)) , (40a) 

PW) = J d q\qa) (-idn(d) 2 r d+1 -^r(q)^j , (40b) 

7 a |V) = Jdq \q a ) ((lT^(q)) ■ (40c) 

Thus q l and 7° act multiplicatively on wave functions, whereas the p, essentially act as 
derivative operators. The commutation relations stated above are clearly satisfied 

A further important point is missing for a successful transition from the classical to the 
quantum theory. The classical covariant Hamiltonian depends on the phase space variables q % 
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and p", so that the Hamiltonian operator would seem to depend on the operators p" for 
which we have not yet given a representation. However, it is always possible to replace these 
operators by pi, as we will now show. Acting on wave functions the pf are realized by 

#~-„(d)j-"V^> (4i) 

which may be derived from an application of our quantization postulates to the classical 
Poisson bracket {q l ,Pj} = t a Sj. It follows that any occurrence of p" can be replaced by 

^-A^'' <42) 
so that the quantum Hamiltonian becomes effectively a function of the operators q l and pY 
On one-dimensional E one reobtains p\ = pi as in the classical theory. 

In the case where the fields q l form the coordinates of a target space manifold M, one must 
take care of appropriate integration measures in the state expansions, and of the fact that 
8(q — q) is a density. The appropriate Schrodinger representation would contain covariant, 
not partial, differentiation operators on M. 

E. Hilbert space and probability interpretation 

The E-dependent states \^(a a )) are elements of a Hilbert space H. The essential algebraic 
structure on a Hilbert space is a scalar product, i.e., a bilinear form (-|-) : H x H — > R. We 
may define such a scalar product in terms of the wave functions corresponding to Hilbert- 
space states: 

= -i J dqi><f) = ~ij dq^r ^. (43) 

Note that the appearance of the gamma matrix T of the local Lorentzian tangent spaces 
guarantees that the scalar product maps to a diffeomorphism scalar of E. The spinor indices 
of the wave functions are suppressed, their E- and field-dependence. 

The necessary requirement that the scalar product should return a scalar function on E 
results in its indefiniteness: indeed, 

(^) = -i j dq^T°ij ^ J dq^ij, (44) 

where the Dirac matrix mixes the spinorial components to prevent a generically positive 
result. The indefiniteness is a feature of quantization on manifolds E of dimension d > 1. 
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For d = 1, the non-equality above becomes an equality (and ip> is simply the complex 
conjugate). 

One of the consequences of this construction is the following: the self-adjoint operators 
with respect to our scalar product are no longer Hermitian with = O. By definition, 
self-adjoint operators satisfy (pip | 0) = (-0 | 00) for all ip and 0. Here this implies O is 
self- adjoint, if, and only if, 

6* = -r°d^r° = 6 . (45) 

But self-adjoint operators are not guaranteed to have real eigenvalues because of the in- 
definiteness of the scalar product, which may produce null states with {ip\ip) = 0. Another 
condition is needed: a self-adjoint operator is orthogonally diagonalizable with real eigenval- 
ues if it is a so-called Pesonen operator satisfying (ip \ dip) 7^ for all null states if), see j^J. 
For the special case d = 1, we have T a = i such that the relation above selects Hermitian 
O, as is the case in conventional quantum mechanics. 

The standard interpretation of quantum mechanics interprets the squared modulus \ip\ 2 
of the Schrodinger wave function ip as a probability density. This is enabled by the fact that 
the Schrodinger equation guarantees the constancy of the total probability j dxip^ip in time. 
To give a similar interpretation here, we need a similar statement. Since the square {ip\ip) is 
no longer positive-definite for d > 1 and hence no longer admits a probability interpretation, 
we need to find another quantity that does. Consider the following vector current on E, 

j a = - dq-p^ip. (46) 



We assume that our covariant Hamiltonian has essentially real eigenvalues, meaning that it 
is self-adjoint with H# = — T°H^T° = H. This is consistent because H is directly related 
to 7 a V a , and we have 7 a ^ = r°7 c T° and Hermitian zV a . Then the Dirac form of relation 
(Jinj) is sufficient to prove that j a is conserved, 

V a j a = 0. (47) 

Locally, a conserved current implies a conserved charge: in normal coordinates on E, 

d a o J d d ~ x a f = . (48) 

Thus we can interpret the integral of j° as the total probability to find the fields in any 
configuration q anywhere on the spatial part of E. This also means that 

p(a a ; q) = -^(a a ; g) 7 (a)^( ( r- q) , (49) 
16 



gives the probability density of finding the field configuration q at a given point with coor- 
dinates a a of S. In normal coordinates, the probability density becomes — ipT ^ = ip'ip, and 
is positive definite. In the one-dimensional case S = M, we can always find global normal 
coordinates, so that we once again recover conventional quantum mechanics. 

The wave function ip(a a ; q) contains all the necessary information to reconstruct the 



Schrodinger wave functional o 
exists, as has been shown in 



the conventional canonical approach in cases where the latter 



25, 



261 ] . This suggests that the covariant canonical approach 



is at least as powerful as the conventional one; more powerful, in fact, since it works on 
backgrounds that do not allow the conventional spacetime foliation by spatial hypersurfaces. 



F. Classical limit and Ehrenfest equations 

To see how the classical limit emerges from our formalism, let us consider the following 
Ehrenfest-type theorem, arising from the expectation value of the Heisenberg equation of 
motion (}32|) above: 

n(d)r d+1 V a <7°/> = <[/,^]>- (50) 

Because the inner product is picture-independent, this equation in particular holds in the 
Schrodinger picture. Assuming now a classical covariant Hamiltonian of the form 

H = --p«pl + V(q l ) (51) 

for constant a, one finds H = ap l pi/ (2dn(d) 2 ) + V(q l ). We wish to evaluate the theorem for 
/ i — ► p> as well as for / i— > 7 a g- J . We use the commutator relations 

BV 

[p„H] = -idnidfr^—iq 1 ), (52a) 
M,H] = ial~ d+l la fP (52b) 



to arrive at the following equations: 

idB_ 



= - v ^ a >> ( 53a ) 
(fl(«*)) = v ^'>- ( 53b ) 

Therefore, the classical covariant Hamiltonian field equations ((7J) are fulfilled within expecta- 
tion values in the form of Ehrenfest equations. Interestingly, no choice of the normalizations 
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of the length scale / or of n(d) was required to obtain this classical result, showing once more 
the consistency of the theory. 



IV. ELEMENTARY APPLICATIONS 



In this section we will give several simple applications of the covariant canonical quantiza- 
tion method. This demonstrates the technique in some detail, but, more importantly, yields 
a number of interesting results: the Klein-Gordon equation arises as the wave equation of 
the relativistic point particle without any need to refer to the relativistic energy-momentum 
relation; the quantization of any bosonic field on an extended parameter space £ with d > 1 
creates spinorial wave functions; in particular, the Dirac equation emerges from the Klein- 
Gordon Lagrangian, along with the prediction of a fermion mass gap and a hierarchy of 
fermions that differ only by their masses. 



A. Relativistic point particles the Klein-Gordon equation 

The simplest application of our quantization formalism is to the relativistic mechanics of 
a point particle. Completely side-stepping its usual derivation from the relativistic energy 
relation, we will find that the quantum wave equation is the Klein-Gordon equation. 

Consider the following action for fields x l : R — > M 13 which describe the worldline 
embedding into a flat Minkowski spacetime, 

S = J doy[=gZ l - (m 1 g aa V <T x i V t ,x i rHj + m 2 ) . (54) 

Variation yields the equation of motion d u {^—g ua g' T,y d a x l ) = 0, and also the gravitational 
constraint g aa V ' a x % V ' a x' J 'rjij = m 2 . Both of these equations are needed to show that the 
above action is classically equivalent to the standard action of the relativistic point particle, 
m J da y' '—d a x i d a x' J 'r/ij , for two mass parameters mi and m 2 which satisfy mim 2 = m 2 . 

The covariant momenta are = m 1 g aa 'V CT x : 'riji, and yield the covariant Hamiltonian 
H = ~2m^9 aa PiPj r l 1 '' ~ \ m 2- Using the Schrodinger representation of the momentum oper- 
ators on wave functions returns the quantum wave equation 

1 _□ _ H^l _ idj) iP(a; x l ) = , (55) 



2mi 2 
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where the box denotes the d'Alembertian on M 1)3 . We still have to deal with the gravitational 
constraint on the classical momenta: in its quantum version it reads (□+m 2 )-?/> = 0. This can 
be satisfied consistently with the wave equation by choosing d a ip{a\ x l ) = 0. The resulting 
wave function is then automatically independent of a, and the equation for ip(x l ) becomes 
the Klein-Gordon equation on Mi j3 . 

The quantum wave equation is a Schrodinger equation with time parameter a, but inde- 
pendence of the wave function of this parameter is forced by the constraint. In this sense, 
classical reparametrization invariance directly implies the Klein-Gordon equation. 

B. Free bosonic strings — Weyl spinors 

One of the most characteristic features of covariant canonical quantization is the fact 
that any bosonic field on a parameter space of dimension d > 1 produces spinor wave 
equations. To illustrate this point in the simplest setup, we consider free bosonic strings on 
a flat target space, given as maps E — > M l n from the two-dimensional worldsheet £ into 
Minkowski space M ljTl . We employ the Polyakov action 

S = ~J d 2 a^ ±g ab d a X%X3 VtJ . (56) 

which is classically equivalent to the Nambu-Goto action that measures the area of the string 
worldsheet, if the gravitational constraint g ab g cd r]* d = r]* cd , where r/* b = d a X l d b X : >ri i j is the 
pull-back of i] to the worldsheet, is implemented. The covariant momenta are derived as 
Pi = —d a Xi and give rise to H = —PiP l a /2. This leads to the Schrodinger picture wave 
equation 

(/- 1 n(2)D + ^da) 4>(a b ; X") = , (57) 

where □ denotes the target space d'Alembert operator on M 1>n .It is illustrative to use the 
separation ansatz ip = ^(a a )^(X^) with spinorial \P for the wave function. This introduces 
a separation constant M, and generates the two separate equations 

(□ + /Mn(2)' 1 ) $ = 0, (58a) 
( 7 a <9 a - M) * = . (58b) 

Thus the quantization of strings generates Weyl spinors on the two-dimensional worldsheet, 
whose mass M is linked to a Klein Gordon equation on the target spacetime. 
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As noted above, it is necessary to satisfy the gravitational constraint. Classically, we may 
rewrite it as p m p\ — g ab p1p l c /2 = which is symmetric. Keeping the symmetry we hence find 
the Schrodinger representation 

{-l {a l b) +9°*)p% = 0- (59) 

The expression in brackets is identically zero by the properties of the Dirac algebra. Surpris- 
ingly the constraint is satisfied automatically in the quantum theory. However, the quanti- 
zation of the string in this formalism has no discernible relation to string theory, where the 
quantum requirement of the gravitational constraint leads to the Virasoro algebra. 

C. Free scalar fields the Dirac equation 

The free scalar field : S 13 -^Rona four-dimensional Lorentzian spacetime S 13 with 
metric g is governed by the Lagrangian 

L = - l -g ab d a <i ) d^- l -m 2 ( f. (60) 

The covariant momenta follow as it a = —d a (f) and yield H = —7i a 7i a /2 + m 2 2 /2. From 
quantization we hence obtain the Schrodinger picture wave equation 

(-^ + |^ 2 -^)^) = - < 61 > 

To illustrate the consequences of our theoretical construction, it is useful to consider a 
separation ansatz tp = \l/(<T a )$(0) for the wave function, where only ^ is taken spinorial. 
This introduces a separation constant M, and the wave equation generates two separate 
equations of the form 

( 7 a <9 a -M)^ = 0. (62b) 

Thus the spacetime dependence of the wave function is described by a Dirac equation with 
mass M. However, this mass is not unconstrained: a spectrum of allowed masses is generated 
by the first equation, which is just the Schrodinger equation of a one-dimensional harmonic 
oscillator with frequency uj = 2m. The mass spectrum is therefore given by 

M k = m (2k + 1) (63) 
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with integer k > 0. While this spectrum does not look particularly appealing phenomeno- 
logically, it should be noted that we obtained the prediction of a mass hierarchy of otherwise 
identical Dirac particles from the non-interacting Klein-Gordon Lagrangian only. With the 
addition of interactions to the scalar action, more complex mass spectra could be generated, 
leading to the possibility of obtaining the generations of the Standard Model from a suitably 
tuned interacting scalar Lagrangian. Another important prediction from equations (|62|) is 
the existence of a mass gap for the fermions: M = is generally not a solution if the 
harmonic potential V(4>) ~ 2 is replaced by a generic potential V((f>). 



D. Local gauge invariance — gauge fields 

We have seen that the quantization of pure scalar field models generates fermionic parti- 
cles with an allowed mass spectrum given by the covariant Hamiltonian of the scalar field. 
Phenomenological relevance additionally requires these fermions to be charged. A mecha- 



nism to serve this purpose has been identified by Weyl [27[ a long time ago, and we will now 
demonstrate its effect. 

The basic observation underlying this mechanism is the invariance of the interpretation- 
ally relevant probability current j a = — i (ip | 7 a | ip) defined in (|4fi|) under local phase shifts 
of the wave function, so-called gauge transformations. These are transformations ip i— > e ip 
under a function A : £ — > R. The quantum evolution equation ()26|). however, is only invari- 
ant under global gauge transformations with constant A. If, in addition, local invariance is 
required, we have to amend this equation by the introduction of a gauge field: 

(H - n{d)l- d+l {f - ieA )) |V>) = . (64) 

It now follows that if is a solution of this equation, then so is the locally gauge transformed 
e ieA ag j on g ag gauge field transforms at the same time as A a \— > A a + d a A. The 
resulting equation is essentially the Dirac equation for particles of mass (H) and charge e. 

To make this statement more precise, observe that in order to speak about the spectrum 
of the theory we do not require knowledge about the probabilities for the original scalar 
fields, but only about the generated fermions. This means we may consider the integrated 
expectation value of the quantum evolution equation as a functional of i/j, 

Sty) = J d d a(ij\ (-n{d)-H d - l H + ¥ - xeA^j |V) (65) 
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in which the wave function's dependence on the original scalar fields is effectively inte- 
grated out. Indeed, if we assume the wave equation has been solved by a product ansatz 
■ip = \l/(cr a ) ( l ) (0) as in the previous examples, we find 

S 9 [V] =(^J d(f)&<^ J d d a^ (V - iefi - M)V , (66) 

for one of the rescaled mass eigenvalues M = n(rf)~ 1 / d_1 J d<fi&H&/ J <i0$f$ in the mass 
spectrum generated by the scalars' covariant Hamiltonian. So integrating out the original 
scalar field freedom, because it is not observable, returns precisely the Lagrangian theory 
for the Dirac field with mass M. 

Now consider a multiplet of N scalar fields fa with a Lagrangian which is invariant under 
the action of some nonabelian subgroup G of SO(N). In this case, the covariant Hamiltonian 
will inherit the G-invariance, resulting in a degeneracy in the spectrum of the theory. The 
energy levels M n can then be labelled by the irreducible representations of G, with the 
degeneracy of each level given by the dimension d n of the irreducible representation under 
which it transforms. A solution of the quantum evolution equation can then be decomposed 
in terms of eigenfunctions of the covariant Hamiltonian as 

oo d n 

1>(<T, <f>) = E *»,«(*)*n,a(0) (67) 

n=0 a=l 

Putting this decomposition into the definition of S[ip] above and using the orthogonality 
relation 

/ d<l&l a (<i>)$ m> f ) (<i>)=6 m>n 6 a ,f ) (68) 

J M 

for the wavefunctions, we arrive at the action for \l/ 

oo „ 

Sm = V / daV n , a (f - M n ) K (69) 

n=0 J Z V 7 

where each kind of fermion has a mass M n and is invariant under a SU(rf n ) symmetry in 
addition to the U(l) symmetry discussed above. Using the same arguments as before, these 
SXJ(d n ) symmetries should also be gauged, giving rise to a spectrum of nonabelian gauge 
symmetries. 

In the case of a finite scalar symmetry group G, there is only a finite number of irreducible 
representations, and correspondingly the gauge group of the fermionic theory then is a finite 
product of SU(iV) gauge groups. For continuous G, an infinite product of SU(iV) factors 
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ensues. Groups that could give a Standard Model-like gauge group SU(3)xSU(2)xU(l) 
(although we haste to point out that this simple model does not give the chiral couplings of 
the Standard Model) include the point groups Td and O. 

Hence, covariant canonical quantization with the additional requirement of local gauge 
invariance is able to produce, from a scalar field Lagrangian, all particles (so far) observed 
in Nature, namely fermions and gauge fields. While their masses are constrained by the 
covariant Hamiltonian of the original classical scalar field theory, the charges are additional 
input at the quantum level of the theory. In order to obtain a more complete theory, 
appropriate gauge-invariant dynamics for the gauge fields have to be added at this stage, 
which is a freedom that we have. 



V. DISCUSSION 



Considering fields as maps from a parameter space X to a target space M, we have 
constructed a covariant quantization method that keeps the diffeomorphism invariance be- 
tween the parameters of £ intact. Covariant canonical quantization is based on the classical 
Hamiltonian theory developed by de Donder jlj and Weyl |2J which makes use of a finite- 
dimensional multi-symplectic phase space, where every field has a set of conjugate momenta 
associated to each of its partial derivatives. The classical theory is completely equivalent 
(in the sense of generating the same solutions) to the conventional Hamiltonian point of 
view, avoiding, however, the explicit breaking of diffeomorphism invariance that arises from 
singling out a time coordinate normal to an assumed foliation of £ by spatial hyper surf aces. 

We have introduced the notion of a covariant Poisson bracket within the classical theory, 
before applying two well-motivated quantization postulates. The first postulate replaces, 
according to Dirac's argument, the Poisson bracket of phase space functions by the com- 
mutator of corresponding operators acting on some Hilbert space. The second postulate 
is geometrically motivated: on parameter spacetimes S of dimension d > 1 it introduces 
the Clifford algebra of Dirac matrices into the quantum theory. The construction is such 
that covariant canonical quantization coincides with conventional canonical quantization in 
d — 1, i.e., when the fields depend only on time. 

The two quantization postulates, applied to the classical de Donder-Weyl theory, for the 
first time allow the derivation of a quantum evolution equation in terms of the covariant 
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Hamiltonian, which had been conjectured before on the mere basis of analogies jal2|- This 
evolution equation effectively unifies the Dirac and the Schrodinger equations. We have 
further developed the theory in the covariant Schrodinger picture, including a discussion 
of the representation of the field and multi-momentum operators, and the relevant Hilbert 
space. Diffeomorphism invariance requires an indefinite inner product on the Hilbert space 
whose consequences for the diagonalizability of operators and their eigenvalues have been 
discussed. We have also provided a probability interpretation for the fields' wave function. 

Further development of the theory could progress in several directions. One of the 
obvious questions concerns the Heisenberg picture for the theory. The formal solution 
for the evolution operator follows from equation ()29j) as the path-ordered exponential 
U(a,(To) ~ V exp (l d ~ l / (dn(d)) j^H^/ad^ ^ . As such it is only locally defined, and path- 
dependent, at least on generic curved spacetimes E where we had to specify the path along 
which the exponent is integrated by using geodesic normal coordinates on sufficiently small 
neighborhoods. This raises questions about the possibility of developing scattering theory 
in these cases. Another question is that of the quantization of gravity. Although we have 
made use of gravitational constraints in two of the examples, it is not so clear how the 
theoretical setup could be consistent with additional dynamics for the background metric 
on the parameter space S. 

We have discussed a number of elementary applications of the formalism with some sur- 
prising and very interesting results. The quantization of the relativistic point particle imme- 
diately yields as quantum wave equation the Klein-Gordon equation, without conceptually 
employing the relativistic energy momentum relation. Though there is no readily discernible 
connection to string theory, it is interesting to note that the quantization of bosonic strings 
produces Weyl spinors on the worldsheet whose mass is linked to a Klein-Gordon equation 
on the target space. This fact expresses one of the most characteristic features of covariant 
canonical quantization, namely that the quantization of any field on £ with d > 1 produces 
a spinorial wave equation. 

We have shown that a purely scalar classical Lagrangian can produce, upon quantization, 
a theory of Dirac fermions interacting with gauge fields. The latter come into the theory 
by requiring the local gauge invariance of the interpretationally relevant probability current 
also on the level of the wave equation. This means that the basic equations of the quantum 
field theory of the Standard Model may emerge from a classical scalar field theory as wave 
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equations. Several interesting results are obtained along with this mechanism: 

The emergence of spinor fields as a purely quantum phenomenon bypasses the usual 
need for a semi-classical treatment of the Dirac equation, which emerges as a fully quantum 
equation from the very beginning. More intriguingly, on the phenomenological side, the 
procedure of covariant canonical quantization provides a new mechanism to unify particles 
of different masses, which leads to an (at least qualitative) explanation of the generations 
of the Standard Model in terms of a self-interacting scalar field. It should be noted that 
in this framework fermionic masses are generated without a Higgs mechanism. Instead, the 
mass and self-interaction of the underlying scalar field manifest themselves by generating a 
mass spectrum for the effective fermion field which is the spacetime part of the scalar field's 
quantum wave function. 

It is interesting to speculate where a scalar model, complicated enough to yield the 
standard model upon quantization, could come from in the first place. One possibility seems 
to be the compactification of a higher- dimensional bosonic, maybe gravitational, theory. 
Such compactifications generally produce a large number of scalar fields as shape moduli 
of the internal manifold, with self-interactions through some effective potential. If such a 
potential had a minimum, one would expect it at some negative value, due to a geometric no- 



go theorem 



28j (see also the discussion in |29j), thus generating a discrete mass spectrum for 



the fermionic fields in the quantum theory which could have phenomenological relevance. 
More speculatively still, the fact that the quantization of a purely scalar classical theory 
necessarily leads to a fermionic quantum theory might be indicating some sort of semi- 
classical supersymmetry at work behind the scenes. 

Covariant canonical quantization clearly does not replace conventional quantum field the- 
ory; rather it adds a first quantization to the usual procedure, which then literally becomes 
a second quantization: first, a scalar model generates spinorial quantum wave equations 
which, after integrating out the unobservable scalar degrees of freedom, become the clas- 
sical equations of motion underlying the Standard Model. Then the quantization of these 
equations proceeds in standard quantum field theoretical fashion. The intriguing result of 
this investigation is the fact that the Standard Model's classical equations of motion for 
fermions and gauge fields, along with the prediction of discrete mass spectra of identically 
charged particles, are generated from the quantization of a purely scalar classical field the- 
ory. This gives reason to hope that covariant canonical quantization might find further uses 
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and applications. 
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